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Abstract. In this letter we present a general classification of integrable models 
of identical classical spins coupled via the isotropic Heisenberg Hamiltonian. Our 
constructive proof of integrability provides a solution scheme for the equations 
describing the time evolution of each spin. We show that for a broad class of 
Heisenberg couplings, it is possible to solve the equations of motion iteratively by 
building a hierarchy of spins, each of which precesses about the composite spin of 
the previous level. This leads to the identification of the frequencies of precession 
at each level. 
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1. Introduction 

The study of the dynamics of classical Heisenberg spins has long proved invaluable 
in a variety of physical contexts. From a mathematical point of view, spin clusters 
have been a natural testing ground to explore the emergence of chaos [1]. The recent 
developments in the experimental study of magnetic molecules [2-4] have generated a 
renewed interest in the dynamics of classical Heisenberg spins. Integral representations 
of the time-dependent spin-spin correlation functions have been presented for a variety 
of small integrable clusters [5-8] . Analogous results for the equivalent neighbor model 
or for an integrable deformation of it have also been obtained [9,10]. Among the more 
recent developments, in [11] the neutron scattering structure factors for a wide range 
of finite quantum spin systems were computed analytically by Haraldsen, Barnes, and 
Musfeldt. 

In this letter we report on a broad classification of integrable cases. In section [21 
we present a class of models leading to integrable equations of motion. In section [3] 
the integrability of the models is proved and the precession frequencies are obtained. 
Section 21 contains our conclusions. 
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2. Classification of the integrable clusters 



The most general Heisenberg Hamiltonian for N classical spins normalized according 
to \Si\ = 1, i = 1,..., N, is 



H 



1 



N 



(Jij Si ■ Sj — S ' gij Si ■ Sj 



(1) 



The matrix of exchange couplings [gij] in is symmetric, g^ = gji, and gn — 0. The 
time evolution of each spin is determined by the differential equations 



N 

9ij Si x Sj 



(2) 



The structure of the equations of motion J2J) explicitly shows that the norm of each 
spin is indeed constant. The equations J2J are non-linear and, for an arbitrary matrix 
of couplings G = [gij], they do not admit a general solution in closed form. In the 
following we identify a class of coupling matrices for which the general solution can 
be explicitly constructed. 

Let us partition the set / = {1,2,3, .. . , N} of spin labels into nonintersecting 
subsets A, B, C, . . . , such that 



AUBUCU. 



AnB 



Anc 



(3) 



We are going to call the set of spins whose labels are in A the 'spin cluster A,' and so on 
for every other subset of /. It is useful to introduce also the following terminology. A 
coupling gij is called internal if its indices i and j belong to the same cluster, whereas 
a coupling gij is called external if i and j belong to different clusters. 

As we shall prove in section [21 if there is a partition J3J) of /, for which all the 
external couplings are equal, the problem of solving the equations of motion (0) is 
reduced to separately solving the spin systems described by the index sets A, B, etc. 
If all the external couplings are equal to J, the coupling matrix G has the form 



G = 



/ G A 
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(4) 



where G A is the matrix of the internal couplings between spins in cluster A, G B the 
matrix of the internal couplings between spins in cluster B, and so on. E is a matrix 
of appropriate dimension, all elements of which are equal to 1. 

Moreover, we claim that if each of the coupling matrices G A , G B , . . . , can in 
turn be written in the form with a single coupling Ja among the sub-clusters of 
A, a single coupling Jb among the sub-clusters of B, etc., (Ja, Jb, ■ ■ ■ , need not be 
equal to J) , and if this can be carried on until we get clusters containing only single 
spins, the equations of motion Q can be solved exactly. As an example, an 8-spin 
model of this type is shown in figure 1. We also notice that all the models studied 
in [5-10] fit into this scheme. 





Figure 1. A model of 8 interacting spins with 4 different couplings. Spins Si 
and S2 form a dimer with a coupling J4. Spin Ss is coupled to this dimer with 
a coupling J3 to form an isosceles triangle. Spins S4, S5, Sa, and S7 form a 
tetrahedron with coupling J2. These individual clusters, together with the single 
spin S3, are coupled to one another by the external coupling J\. 



3. Solution method and precession frequencies 

If the coupling matrix has the form ijlf. the equation of motion for a spin belonging 
to cluster A is 

£ 9aa> S <*xS a ' + JS a xS A , VaeA (5) 

a'eA 

where S A = Si is the sum of all the spins outside cluster A. Similar equations 

iel\A 

can be written for the spins in every other cluster B, C, etc. 
With the definitions 

Sa = X , Sb = X Sp, ... (6) 

a£A P£B 

we obtain the equations of motion describing the time evolution of the spin sums © 
by adding the equations |5j) for all spins in a given cluster. E.g., for cluster A we 
obtain 

S A = JS A x S A . (7) 
Since the total spin S of the system can be written as 

S = S A + S A (8) 
equation becomes 

S A = JS A x S . (9) 
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Since the total spin S is conserved in the Heisenberg model, equation 10 can be 
solved immediately. The solution describes the precession of the vector Sa around S. 
Choosing a coordinate system such that S = Sz, the solution can be written as 

S A (t) = S l [z + Si [xcos(wt - cj> A ) - y$m{ut - (j> A )\ (10) 

where lu = JS, S A is the projection of S A along the direction of S, S A is its projection 
perpendicular to the direction of S, and 4> A is the initial phase of S A . Since all the 
S A (t) and S A (t) = S — S A (t) are now explicitly known functions of time, solving 
equation JSJ is reduced to solving the internal coupling sector for each cluster. Indeed, 
let us focus again on cluster A and let us assume that its index set A can itself be 
partitioned in a way similar to J2J), i.e., 

A = MUNUPU... , MniV = 0, MHP = 0, (11) 

and that the coupling matrix G A has the structure 10} , with all external couplings 
between the sub-clusters M, N, ... , being equal to J A . The equations of motion 
involving the spins of sub-cluster M can then be written as 

S„ = J2 x + J aS„ x S M + JS„ x S A , V/i e M (12) 

where S M = Si. Once again, equations of the same form will hold for the spins 

i£A\M 

in the other sub-clusters of A. 

Adding all the equations for the spins in sub-cluster M in a way similar to the 
derivations of (Q) and ©, and defining Sm = J]] S^, we obtain 

fj,eM 

Sm = JaSm x S A + JSm x S a (13) 

which can also be written in the form 

S M = JS M x S + (J A - J)S M x S A . (14) 

As long as we can apply the reduction scheme of Section 2 to the next level of 
sub-clusters, this solving procedure can be carried on iteratively. For the purpose of 
illustration, suppose that the internal couplings between the spins in sub-cluster M are 
all equal to Jm, i-e., that the reduction procedure ends at this level.f The equations 
of motion for the individual spins in cluster M can then be written as 

Si = JSi x S + {J A - J)S t x S A + {Jm - Ja)S 1 x S m , i e M. (15) 

In order to write the solution for Si(t) in a closed form we apply successively the 
Fourier transform method used in [6,7,10] to find S A (t), then Sm(^), and finally Si(t), 
i G M. Though cumbersome and tedious, this is straightforward and mechanical. For 
the purposes of this letter, we restrict ourselves to stating the resulting frequencies of 
the various spin precessions : JS, (J A — J)S A , and (Jm — Ja)Sm- Figure 2 illustrates 
the precessions of the individual spins and the spin composites for the 8-spin model 
in figure 1. 

More generally, if the reduction scheme leads to a hierarchy of n levels of sub- 
clusters, i G c c . . . C 1^ C 1^ C /, where each fc-level cluster has a 

X This is the same as saying that the sub-clusters of M consist of individual spins, each coupled to all 
the others with coupling Jm , or that M is an equivalent neighbor model with respect to its internal 
couplings. 




Figure 2. The time evolution of the 8-spin model of figure 1. The net spin of the 
isosceles triangle, Si28i the net spin of the tetrahedron, S4567, and the individual 
spin S3 precess with angular speed JiS about the total spin of the system, S. 
The individual spins Si, S5, Se, and S7 in turns precess about S4567. The dimer 
spin S12 and the individual spin S$ precess about the isosceles spin Si28- Finally, 
Si and S2 precess about S12. 

coupling J(£) among its (k + l)-level sub-clusters, the time evolution of the i th spin 
St can be found by iteratively solving the equations 

s =0 

S(i) = JS(i) x s 

S(2) = JS (2 ) xS + (J (1) - J)5 (2) x S {1) (16) 

S(n) = JS(n) X S + (J(!) - J)5(„) X 5(i) + . . . + (J( n -1) ~ <^(n-2))5( n ) x ^(n-l) 
5j = JSi X S + (J(!) — J)5, X 5(i) + . . . + (J( n ) - J( n _i))5j X 5(„) 

The Fourier transformed version of the hierarchy of equations (|16(l is: 
iwS = 
iuS^ = JS(i^ (g> 5 

iwS (2) = JS(2) ®S + ( J ( i) - J) 5(2) ® 5 ( i) (17) 

«w5( n ) = J5(„) <g> 5 + (J (1 ) - J)5(„) (8) 5(i) + . . . + (J( n -1) ~ J(n-2))S(n) ® ^(n-1) 

iw5i =JSi®S+ (J ( i) - J)5i ® 5(i) + . . . + (J (n) - J (n -i))5 i 5 (n) 
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where the symbol <£> denotes a vector product in the real space and a convolution in 
the Fourier space. 

Although the explicit solution to the general case (|17f) would be too cumbersome 
to write in a closed form, it is easy to see that the first equation of l|17[) will lead to 
S(lo) oc S(lu) (i.e., o->( ) = or S 1 is constant). The second equation would introduce 
a new frequency in the Fourier spectrum of S^y. u>(i) = JS (for instance, this 
gives the spin precession frequency for each of spins in a dimer [5,8]). The third 
equation will give rise to the additional frequency lu( 2 ) — (J(i) ~ J)S(i) m the Fourier 
spectrum of (as it happens for the individual spins on the sites of a squashed 
equivalent neighbor model [7, 10].) We are therefore led to state that the solution for 
the individual spin Si will contain (^-functions of various sums or differences of the 
frequencies JS, (J(i) - J)S(i), (J( 2 ) - J(i))S( 2 ), ■ ■ ■ , (J( n ) ~ J(n-i))S( n ) or, in other 
words, the Fourier spectrum of Si(t) will contain only these frequencies. 

4. Conclusions 

In this letter we have identified a class of integrable spin cluster models and have 
presented a procedure for obtaining the exact solutions for the time evolution of the 
individual spins, listing the frequencies that appear in their Fourier spectra. We hope 
that these results and solving procedures will prove useful to anyone interested in 
easily identifying and solving specific classical Heisenberg spin models relevant for 
describing the physical properties of magnetic molecules. 
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